Preclinical exploration of the DNA damage response pathway using the interactive neuroblastoma cell line explorer CLEAN.

Neuroblastoma (NB) is the most common cancer in infancy with an urgent need for more efficient targeted therapies. The development of novel (combinatorial) treatment strategies relies on extensive explorations of signaling perturbations in neuroblastoma cell lines, using RNA-Seq or other high throughput technologies (e.g. phosphoproteomics). This typically requires dedicated bioinformatics support, which is not always available. Additionally, while data from published studies are highly valuable and raw data (e.g. fastq files) are nowadays released in public repositories, data processing is time-consuming and again difficult without bioinformatics support. To facilitate NB research, more user-friendly and immediately accessible platforms are needed to explore newly generated as well as existing high throughput data. To make this possible, we developed an interactive data centralization and visualization web application, called CLEAN (the Cell Line Explorer web Application of Neuroblastoma data; https://ccgg.ugent.be/shiny/clean/). By focusing on the regulation of the DNA damage response, a therapeutic target of major interest in neuroblastoma, we demonstrate how CLEAN can be used to gain novel mechanistic insights and identify putative drug targets in neuroblastoma.

Benchmark of tools for in silico prediction of MHC class I and class II genotypes from NGS data.

BACKGROUND: The Human Leukocyte Antigen (HLA) genes are a group of highly polymorphic genes that are located in the Major Histocompatibility Complex (MHC) region on chromosome 6. The HLA genotype affects the presentability of tumour antigens to the immune system. While knowledge of these genotypes is of utmost importance to study differences in immune responses between cancer patients, gold standard, PCR-derived genotypes are rarely available in large Next Generation Sequencing (NGS) datasets. Therefore, a variety of methods for in silico NGS-based HLA genotyping have been developed, bypassing the need to determine these genotypes with separate experiments. However, there is currently no consensus on the best performing tool. RESULTS: We evaluated 13 MHC class I and/or class II HLA callers that are currently available for free academic use and run on either Whole Exome Sequencing (WES) or RNA sequencing data. Computational resource requirements were highly variable between these tools. Three orthogonal approaches were used to evaluate the accuracy on several large publicly available datasets: a direct benchmark using PCR-derived gold standard HLA calls, a correlation analysis with population-based allele frequencies and an analysis of the concordance between the different tools. The highest MHC-I calling accuracies were found for Optitype (98.0%) and arcasHLA (99.4%) on WES and RNA sequencing data respectively, while for MHC-II HLA-HD was the most accurate tool for both data types (96.2% and 99.4% on WES and RNA data respectively). CONCLUSION: The optimal strategy for HLA genotyping from NGS data depends on the availability of either WES or RNA data, the size of the dataset and the available computational resources. If sufficient resources are available, we recommend Optitype and HLA-HD for MHC-I and MHC-II genotype calling respectively.

ALK fusion NSCLC oncogenes promote survival and inhibit NK cell responses via SERPINB4 expression.

Anaplastic lymphoma kinase (ALK) fusion variants in Non-Small Cell Lung Cancer (NSCLC) consist of numerous dimerizing fusion partners. Retrospective investigations suggest that treatment benefit in response to ALK tyrosine kinase inhibitors (TKIs) differs dependent on the fusion variant present in the patient tumor. Therefore, understanding the oncogenic signaling networks driven by different ALK fusion variants is important. To do this, we developed controlled inducible cell models expressing either Echinoderm Microtubule Associated Protein Like 4 (EML4)-ALK-V1, EML4-ALK-V3, Kinesin Family Member 5B (KIF5B)-ALK, or TRK-fused gene (TFG)-ALK and investigated their transcriptomic and proteomic responses to ALK activity modulation together with patient-derived ALK-positive NSCLC cell lines. This allowed identification of both common and isoform-specific responses downstream of these four ALK fusions. An inflammatory signature that included upregulation of the Serpin B4 serine protease inhibitor was observed in both ALK fusion inducible and patient-derived cells. We show that Signal transducer and activator of transcription 3 (STAT3), Nuclear Factor Kappa B (NF-κB) and Activator protein 1 (AP1) are major transcriptional regulators of SERPINB4 downstream of ALK fusions. Upregulation of SERPINB4 promotes survival and inhibits natural killer cell-mediated cytotoxicity, which has potential for therapeutic impact targeting the immune response together with ALK TKIs in NSCLC.

A clinically annotated post-mortem approach to study multi-organ somatic mutational clonality in normal tissues.

Recent research on normal human tissues identified omnipresent clones of cells, driven by somatic mutations known to be responsible for carcinogenesis (e.g., in TP53 or NOTCH1). These new insights are fundamentally changing current tumor evolution models, with broad oncological implications. Most studies are based on surgical remnant tissues, which are not available for many organs and rarely in a pan-organ setting (multiple organs from the same individual). Here, we describe an approach based on clinically annotated post-mortem tissues, derived from whole-body donors that are routinely used for educational purposes at human anatomy units. We validated this post-mortem approach using UV-exposed and unexposed epidermal skin tissues and confirm the presence of positively selected NOTCH1/2-, TP53- and FAT1-driven clones. No selection signals were detected in a set of immune genes or housekeeping genes. Additionally, we provide the first evidence for smoking-induced clonal changes in oral epithelia, likely underlying the origin of head and neck carcinogenesis. In conclusion, the whole-body donor-based approach provides a nearly unlimited healthy tissue resource to study mutational clonality and gain fundamental mutagenic insights in the presumed earliest stages of tumor evolution.

ATR inhibition enables complete tumour regression in ALK-driven NB mouse models.

High-risk neuroblastoma (NB) often involves MYCN amplification as well as mutations in ALK. Currently, high-risk NB presents significant clinical challenges, and additional therapeutic options are needed. Oncogenes like MYCN and ALK result in increased replication stress in cancer cells, offering therapeutically exploitable options. We have pursued phosphoproteomic analyses highlighting ATR activity in ALK-driven NB cells, identifying the BAY1895344 ATR inhibitor as a potent inhibitor of NB cell growth and proliferation. Using RNA-Seq, proteomics and phosphoproteomics we characterize NB cell and tumour responses to ATR inhibition, identifying key components of the DNA damage response as ATR targets in NB cells. ATR inhibition also produces robust responses in mouse models. Remarkably, a 2-week combined ATR/ALK inhibition protocol leads to complete tumor regression in two independent genetically modified mouse NB models. These results suggest that NB patients, particularly in high-risk groups with oncogene-induced replication stress, may benefit from ATR inhibition as therapeutic intervention.

Network-Based Analysis to Identify Drivers of Metastatic Prostate Cancer Using GoNetic.

Most known driver genes of metastatic prostate cancer are frequently mutated. To dig into the long tail of rarely mutated drivers, we performed network-based driver identification on the Hartwig Medical Foundation metastatic prostate cancer data set (HMF cohort). Hereto, we developed GoNetic, a method based on probabilistic pathfinding, to identify recurrently mutated subnetworks. In contrast to most state-of-the-art network-based methods, GoNetic can leverage sample-specific mutational information and the weights of the underlying prior network. When applied to the HMF cohort, GoNetic successfully recovered known primary and metastatic drivers of prostate cancer that are frequently mutated in the HMF cohort (TP53, RB1, and CTNNB1). In addition, the identified subnetworks contain frequently mutated genes, reflect processes related to metastatic prostate cancer, and contain rarely mutated driver candidates. To further validate these rarely mutated genes, we assessed whether the identified genes were more mutated in metastatic than in primary samples using an independent cohort. Then we evaluated their association with tumor evolution and with the lymph node status of the patients. This resulted in forwarding several novel putative driver genes for metastatic prostate cancer, some of which might be prognostic for disease evolution.

Loss of RET Promotes Mesenchymal Identity in Neuroblastoma Cells.

Aberrant activation of anaplastic lymphoma kinase (ALK) drives neuroblastoma (NB). Previous work identified the RET receptor tyrosine kinase (RTK) as a downstream target of ALK activity in NB models. We show here that ALK activation in response to ALKAL2 ligand results in the rapid phosphorylation of RET in NB cells, providing additional insight into the contribution of RET to the ALK-driven gene signature in NB. To further address the role of RET in NB, RET knockout (KO) SK-N-AS cells were generated by CRISPR/Cas9 genome engineering. Gene expression analysis of RET KO NB cells identified a reprogramming of NB cells to a mesenchymal (MES) phenotype that was characterized by increased migration and upregulation of the AXL and MNNG HOS transforming gene (MET) RTKs, as well as integrins and extracellular matrix components. Strikingly, the upregulation of AXL in the absence of RET reflects the development timeline observed in the neural crest as progenitor cells undergo differentiation during embryonic development. Together, these findings suggest that a MES phenotype is promoted in mesenchymal NB cells in the absence of RET, reflective of a less differentiated developmental status.

Low immunogenicity of common cancer hot spot mutations resulting in false immunogenic selection signals.

Cancer is driven by somatic mutations that result in a cellular fitness advantage. This selective advantage is expected to be counterbalanced by the immune system when these driver mutations simultaneously lead to the generation of neoantigens, novel peptides that are presented at the cancer cell membrane via HLA molecules from the MHC complex. The presentability of these peptides is determined by a patient's MHC genotype and it has been suggested that this results in MHC genotype-specific restrictions of the oncogenic mutational landscape. Here, we generated a set of virtual patients, each with an identical and prototypical MHC genotype, and show that the earlier reported HLA affinity differences between observed and unobserved mutations are unrelated to MHC genotype variation. We demonstrate how these differences are secondary to high frequencies of 13 hot spot driver mutations in 6 different genes. Several oncogenic mechanisms were identified that lower the peptides' HLA affinity, including phospho-mimicking substitutions in BRAF, destabilizing tyrosine mutations in TP53 and glycine-rich mutational contexts in the GTP-binding KRAS domain. In line with our earlier findings, our results emphasize that HLA affinity predictions are easily misinterpreted when studying immunogenic selection processes.

ALK ligand ALKAL2 potentiates MYCN-driven neuroblastoma in the absence of ALK mutation.

High-risk neuroblastoma (NB) is responsible for a disproportionate number of childhood deaths due to cancer. One indicator of high-risk NB is amplification of the neural MYC (MYCN) oncogene, which is currently therapeutically intractable. Identification of anaplastic lymphoma kinase (ALK) as an NB oncogene raised the possibility of using ALK tyrosine kinase inhibitors (TKIs) in treatment of patients with activating ALK mutations. 8-10% of primary NB patients are ALK-positive, a figure that increases in the relapsed population. ALK is activated by the ALKAL2 ligand located on chromosome 2p, along with ALK and MYCN, in the "2p-gain" region associated with NB. Dysregulation of ALK ligand in NB has not been addressed, although one of the first oncogenes described was v-sis that shares > 90% homology with PDGF. Therefore, we tested whether ALKAL2 ligand could potentiate NB progression in the absence of ALK mutation. We show that ALKAL2 overexpression in mice drives ALK TKI-sensitive NB in the absence of ALK mutation, suggesting that additional NB patients, such as those exhibiting 2p-gain, may benefit from ALK TKI-based therapeutic intervention.

11q Deletion or ALK Activity Curbs DLG2 Expression to Maintain an Undifferentiated State in Neuroblastoma.

High-risk neuroblastomas typically display an undifferentiated or poorly differentiated morphology. It is therefore vital to understand molecular mechanisms that block the differentiation process. We identify an important role for oncogenic ALK-ERK1/2-SP1 signaling in the maintenance of undifferentiated neural crest-derived progenitors through the repression of DLG2, a candidate tumor suppressor gene in neuroblastoma. DLG2 is expressed in the murine "bridge signature" that represents the transcriptional transition state when neural crest cells or Schwann cell precursors differentiate to chromaffin cells of the adrenal gland. We show that the restoration of DLG2 expression spontaneously drives neuroblastoma cell differentiation, highlighting the importance of DLG2 in this process. These findings are supported by genetic analyses of high-risk 11q deletion neuroblastomas, which identified genetic lesions in the DLG2 gene. Our data also suggest that further exploration of other bridge genes may help elucidate the mechanisms underlying the differentiation of NC-derived progenitors and their contribution to neuroblastomas.

Lack of detectable neoantigen depletion signals in the untreated cancer genome.

Somatic mutations can result in the formation of neoantigens, immunogenic peptides that are presented on the tumor cell surface by HLA molecules. These mutations are expected to be under negative selection pressure, but the extent of the resulting neoantigen depletion remains unclear. On the basis of HLA affinity predictions, we annotated the human genome for its translatability to HLA binding peptides and screened for reduced single nucleotide substitution rates in large genomic data sets from untreated cancers. Apparent neoantigen depletion signals become negligible when taking into consideration trinucleotide-based mutational signatures, owing to lack of power or to efficient immune evasion mechanisms that are active early during tumor evolution.

Targeting Filamin A Reduces Macrophage Activity and Atherosclerosis.

BACKGROUND: The actin-binding protein FLNA (filamin A) regulates signal transduction important for cell locomotion, but the role of macrophage-specific FLNA during atherogenesis has not been explored. METHODS: We analyzed FLNA expression in human carotid atherosclerotic plaques by immunofluorescence. We also produced mice with Flna-deficient macrophages by breeding conditional Flna-knockout mice ( Flna o/fl) with mice expressing Cre from the macrophage-specific lysosome M promoter ( LC). Atherosclerosis in vivo was studied by transplanting bone marrow from male Flna o/fl/ LC mice to atherogenic low-density lipoprotein receptor-deficient ( Ldlr-/-) mice; and by infecting Flna o/fl and Flna o/fl/ LC mice with AdPCSK9 (adenoviral vector overexpressing proprotein convertase subtilisin/kexin type 9). Furthermore, C57BL/6 mice were infected with AdPCSK9 and then treated with the calpain inhibitor calpeptin to inhibit FLNA cleavage. RESULTS: We found that macrophage FLNA expression was higher in advanced than in intermediate human atherosclerotic plaques. Flna o/fl/ LC macrophages proliferated and migrated less than controls; expressed lower levels of phosphorylated AKT and ERK1/2; exhibited reduced foam cell formation and lipid uptake; and excreted more lipids. The deficiency of Flna in macrophages markedly reduced the size of aortic atherosclerotic plaques in both Ldlr-/-BMT: Flnao/fl/LC and AdPCSK9-infected Flna o/fl/ LC mice. Intima/media ratios and numbers of CD68-positive macrophages in atherosclerotic plaques were lower in Flna-deficient mice than in control mice. Moreover, we found that STAT3 interacts with a calpain-cleaved carboxyl-terminal fragment of FLNA. Inhibiting calpain-mediated FLNA cleavage with calpeptin in macrophages reduced nuclear levels of phosphorylated STAT3, interleukin 6 secretion, foam cell formation, and lipid uptake. Finally, calpeptin treatment reduced the size of atherosclerotic plaques in C57BL/6 mice infected with AdPCSK9. CONCLUSIONS: Genetic inactivation of Flna and chemical inhibition of calpain-dependent cleavage of FLNA impaired macrophage signaling and function, and reduced atherosclerosis in mice, suggesting that drugs targeting FLNA may be useful in the treatment of atherosclerosis.

Elevated pyrimidine dimer formation at distinct genomic bases underlies promoter mutation hotspots in UV-exposed cancers.

Sequencing of whole cancer genomes has revealed an abundance of recurrent mutations in gene-regulatory promoter regions, in particular in melanoma where strong mutation hotspots are observed adjacent to ETS-family transcription factor (TF) binding sites. While sometimes interpreted as functional driver events, these mutations are commonly believed to be due to locally inhibited DNA repair. Here, we first show that low-dose UV light induces mutations preferably at a known ETS promoter hotspot in cultured cells even in the absence of global or transcription-coupled nucleotide excision repair (NER). Further, by genome-wide mapping of cyclobutane pyrimidine dimers (CPDs) shortly after UV exposure and thus before DNA repair, we find that ETS-related mutation hotspots exhibit strong increases in CPD formation efficacy in a manner consistent with tumor mutation data at the single-base level. Analysis of a large whole genome cohort illustrates the widespread contribution of this effect to recurrent mutations in melanoma. While inhibited NER underlies a general increase in somatic mutation burden in regulatory elements including ETS sites, our data supports that elevated DNA damage formation at specific genomic bases is at the core of the prominent promoter mutation hotspots seen in skin cancers, thus explaining a key phenomenon in whole-genome cancer analyses.

Phosphoproteome and gene expression profiling of ALK inhibition in neuroblastoma cell lines reveals conserved oncogenic pathways.

Anaplastic lymphoma kinase (ALK) is a tyrosine kinase receptor that is a clinical target of major interest in cancer. Mutations and rearrangements in ALK trigger the activation of the encoded receptor and its downstream signaling pathways. ALK mutations have been identified in both familial and sporadic neuroblastoma cases as well as in 30 to 40% of relapses, which makes ALK a bona fide target in neuroblastoma therapy. Tyrosine kinase inhibitors (TKIs) that target ALK are currently in clinical use for the treatment of patients with ALK-positive non-small cell lung cancer. However, monotherapy with the ALK inhibitor crizotinib has been less encouraging in neuroblastoma patients with ALK alterations, raising the question of whether combinatorial therapy would be more effective. In this study, we established both phosphoproteomic and gene expression profiles of ALK activity in neuroblastoma cells exposed to first- and third-generation ALK TKIs, to identify the underlying molecular mechanisms and identify relevant biomarkers, signaling networks, and new therapeutic targets. This analysis has unveiled various important leads for novel combinatorial treatment strategies for patients with neuroblastoma and an increased understanding of ALK signaling involved in this disease.

An antisense RNA capable of modulating the expression of the tumor suppressor microRNA-34a.

The microRNA-34a is a well-studied tumor suppressor microRNA (miRNA) and a direct downstream target of TP53 with roles in several pathways associated with oncogenesis, such as proliferation, cellular growth, and differentiation. Due to its broad tumor suppressive activity, it is not surprising that miR34a expression is altered in a wide variety of solid tumors and hematological malignancies. However, the mechanisms by which miR34a is regulated in these cancers is largely unknown. In this study, we find that a long noncoding RNA transcribed antisense to the miR34a host gene, is critical for miR34a expression and mediation of its cellular functions in multiple types of human cancer. We name this long noncoding RNA lncTAM34a, and characterize its ability to facilitate miR34a expression under different types of cellular stress in both TP53-deficient and wild-type settings.

Clinical response of the novel activating ALK-I1171T mutation in neuroblastoma to the ALK inhibitor ceritinib.

Tumors with anaplastic lymphoma kinase (ALK) fusion rearrangements, including non-small-cell lung cancer and anaplastic large cell lymphoma, are highly sensitive to ALK tyrosine kinase inhibitors (TKIs), underscoring the notion that such cancers are addicted to ALK activity. Although mutations in ALK are heavily implicated in childhood neuroblastoma, response to the ALK TKI crizotinib has been disappointing. Embryonal tumors in patients with DNA repair defects such as Fanconi anemia (FA) often have a poor prognosis, because of lack of therapeutic options. Here we report a child with underlying FA and ALK mutant high-risk neuroblastoma responding strongly to precision therapy with the ALK TKI ceritinib. Conventional chemotherapy treatment caused severe, life-threatening toxicity. Genomic analysis of the initial biopsy identified germline FANCA mutations as well as a novel ALK-I1171T variant. ALK-I1171T generates a potent gain-of-function mutant, as measured in PC12 cell neurite outgrowth and NIH3T3 transformation. Pharmacological inhibition profiling of ALK-I1171T in response to various ALK TKIs identified an 11-fold improved inhibition of ALK-I1171T with ceritinib when compared with crizotinib. Immunoaffinity-coupled LC-MS/MS phosphoproteomics analysis indicated a decrease in ALK signaling in response to ceritinib. Ceritinib was therefore selected for treatment in this child. Monotherapy with ceritinib was well tolerated and resulted in normalized catecholamine markers and tumor shrinkage. After 7.5 mo treatment, the residual primary tumor shrunk, was surgically removed, and exhibited hallmarks of differentiation together with reduced Ki67 levels. Clinical follow-up after 21 mo treatment revealed complete clinical remission including all metastatic sites. Therefore, ceritinib presents a viable therapeutic option for ALK-positive neuroblastoma.